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ABSTRACT 

Precursor emission has been observed in a non-negligible fraction of gamma-ray bursts. The time gap 
between the precursor and the main burst extends in some case up to hundreds of seconds, such as in 
GRB041219A, GRB050820A and GRB060124. Both the origin of the precursor and the large value of 
the time gap are controversial. Here we investigate the maximum possible time gaps arising from the jet 
propagation inside the progenitor star, in models which assume that the precursor is produced by the jet 
bow shock or the cocoon breaking out of the progenitor. Due to the pressure drop ahead of the jet head 
after it reaches the stellar surface, a rarefaction wave propagates back into the jet at the sound speed, 
which re-accelerates the jet to a relativistic velocity and therefore limits the gap period to within about 
ten seconds. This scenario therefore cannot explain gaps which are hundreds of seconds long. Instead, 
we ascribe such long time gaps to the behavior of the central engine, and suggest a fallback collapsar 
scenario for these bursts. In this scenario, the precursor is produced by a weak jet formed during the 
initial core collapse, possibly related to MHD processes associated with a short-lived proto-neutron star, 
while the main burst is produced by a stronger jet fed by fallback accretion onto the black hole resulting 
from the collapse of the neutron star. We have examined the propagation times of the weak precursor jet 
through the stellar progenitor. We find that the initial weak jet can break out of the progenitor in a time 
less than ten seconds (a typical precursor duration) provided that it has a moderately high relativistic 
Lorentz factor T ;> 10. The longer (~ 100 s) time gap following this is determined, in this scenario, 
by the fall-back timescale, which is at the same time long enough for the exit channel to close after 
the precursor activity ceases, thus allowing for the collimation by the cocoon pressure of the subsequent 
main jet, as required. 

Subject headings: gamma rays: bursts — radiation mechanisms: non-thermal 



1. INTRODUCTION 

Weak, softer precursor emission, occurring from several 
seconds to hundreds of seconds before the main gamma-ray 
pulse, is occasionally detected in long gamma-ray bursts 
(GRB). Often there is almost no emission detectable in 
the time gap between the precursor and the main event. 
The earliest detections of precursor events were made with 
the Ginga satellite (GRB900126, Murakami et al. 1991). 
Since then, precursors have been found in many bursts, 
e.g. by BeppoSAX (e.g. GRB011121, Piro et al. 2005), 
HETE2 (e.g. GRB030329, Vanderspek et al. 2004), IN- 
TEGRAL (e.g. GRB041219A, Mereghetti et al. 2004) and 
Swift (e.g. GRB050820A, Cenko et al. 2006; GRB060124, 
Romano et al. 2006; GRB061121, Page et al. 2007). Laz- 
zati (2005) analyzed a sample of bright, long BATSE burst 
light curves in the 200 s period before the detection of the 
GRB prompt emission, and found that a sizable fraction 
(~ 20%) of them show evidence for gamma-ray emission 
above the background between 10 to ~ 200 s (the limit of 
the investigated period) before the main burst. The pre- 
cursors are typically characterized by non-thermal spectra 
which tend to be softer than the main burst, and con- 
tain from 0.1% to a few % of the total event counts. For 
some bursts with known redshifts, the isotropic energies 



of the precursors are determined, and can be as large as 
E iso ~5x 10 51 erg for GRB060124 (Romano et al. 2006; 
Misra et al. 2007). The time gap between the precur- 
sor and the main event ranges from several seconds to a 
few hundreds of seconds (e.g. 200 s for GRB041219A and 
GRB050820A, and 570 s for GRB060124 in the observer 
frame) . 

Theoretical models for precursors have been proposed, 
relating it to the breakout of the main GRB jet. As sug- 
gested by Ramirez-Ruiz et al. (2002a) and Waxman & 
Meszaros (2003), when the jet is making its way out of 
the stellar mantle, a bow shock runs ahead and a strong 
thermal precursor is produced as this shock breaks out 4 . 
In this model, as usual, the main gamma-ray burst is pro- 
duced when the internal shock occurs at some radius i? 7 , 
the delay time between precursor and the main burst be- 
ing At ~ R 7 /2T 2 c, where T is the Lorentz factor of the jet. 
Since the bursts with a precursor also display significant 
;> MeV emission, similarly to other typical long bursts, 
their bulk Lorentz factors must be likewise T ;> 100 (e.g. 
Baring & Harding 1997; Lithwick & Sari 2001), hence the 
typical delay At is less then a second for the typical inter- 



nal shock radii i? 7 ~ 10 — 10 cm (see Meszaros (2006) 
and Zhang (2007) for a recent review). Since the variabil- 
ity timescale is also i? 7 /2r 2 c, this scenario would imply 
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a gap timescale comparable to the variability timescale, 
which is clearly not the case. Morsony et al. (2007) inves- 
tigated the jet propagation through the star with an elab- 
orate numerical simulation and identified three distinct 
phases during the evolution, i.e. the precursor phase, the 
shocked jet phase and the unshocked jet phase. They find 
that a high-pressure cocoon, formed as the sub-relativistic 
jet head makes its way out of the star, is released when 
the head of the jet reaches the stellar surface (which can 
produce a precursor, Ramirez- Ruiz et al. 2002b). Within 
the framework of this model, an observer which is located 
at small or modest viewing angles relative to the jet axis 
would be expected to see first a relatively bright precur- 
sor, then a gap phase during which there is little emission, 
and finally a bright GRB phase when the unshocked jet 
reaches the radius of the stellar surface. 

In §2, we consider first a precursor model such as dis- 
cussed above, which is based on the time delays associ- 
ated with the same jet giving rise both to the precursor 
and to the main burst event. We show that the shocked 
jet phase could not last much longer than ~ 10 s, due to 
the fact that once the jet head reaches the stellar surface, 
a rarefaction wave will form which propagates back into 
the shocked jet component, decreasing the pressure in the 
shocked jet significantly, and as a result, the jet attains its 
ultimate relativistic velocity. Hence, the propagation time 
taken by the unshocked jet to reach the stellar surface will 
be limited to ~ R*/c, which is only from a few seconds 
to ten seconds for Wolf-Rayet progenitors of long GRBs, 
where i?+ is the stellar radius of the progenitor. Thus it 
appears that this model has difficulties in explaining the 
<~ 100 s long gaps seen in some GRBs. In §3 we then sug- 
gest that gaps hundreds of seconds long could arise in a 
fallback collapsar scenario, or in a "type II collapsar" sce- 
nario (MacFadyen et al. 2001), the gap time being related 
to the feeding time of the fallback disk, or the instability 
timescale in the accretion disk. We argue that in such sce- 
narios an initial weak jet, which is still able to break out 
through the star, produces the precursor, while the main 
burst is produced by a subsequent more energetic jet, fed 
by accretion of the fallback gas. The initial precursor jet, 
whatever its formation mechanism, must be weak enough 
not to disrupt the star, so that a fallback disk can still 
form later on. This constrains the luminosity of the pre- 
cursor jet to be lower than E^/tj, where ~ 10 51 erg 
is the binding energy of the star envelope and tj is the 
precursor duration. Motivated by the findings of Woosley 
& Zhang (2007) that a weaker jet propagates more slowly 
inside the jet, in §4 we study analytically the dependence 
of the propagation time of the jet on other jet parameters. 
We confirm here the findings of these simulations, and ob- 
tain some constraints on the precursor jet luminosity and 
Lorentz factors. Some possible formation mechanisms for 
the precursor jet are considered in §5. A summary of the 
results is given in §6. 

2. JET DYNAMICS INSIDE THE STAR AND AT THE 
BREAKOUT 

In the collapsar model (Woosley 1993; Paczyhski 1998; 
MacFadyen & Woosley 1999) , GRBs are caused by rela- 
tivistic jets expelled along the rotation axis of a collapsing 
stellar core. The relativistic jets are due to a black hole 
or a neutron star accretion disk, after the iron core of the 



massive star progenitor has collapsed. 

As the jet advances through the star, it drives a bow 
shock ahead of itself into the star, while the ram pressure 
of the shocked gas ahead of the jet drives a reverse shock 
into the head of the jet, slowing the jet head down to a sub- 
relativistic velocity. Thus there are three distinct regions: 
1) In front of the contact discontinuity between the jet and 
the stellar gas, there is a thin layer of shocked stellar gas 
moving ahead with a sub-relativistic velocity ~ Vh into the 
star; 2) behind the contact discontinuity, there is a shocked 
jet region, where the relativistic jet with Tj 3> 1, is slowed 
down to a velocity <~ Vh by the reverse shock, and 3) and 
below this is the unshocked jet, whose bulk Lorentz factor 
behaves as if it were a free jet. In this unshocked jet, at 
lower radii, the gradual conversion of its internal energy 
into kinetic energy results in Tj = To(r0j/ro6o), until a 
saturation radius at which it reaches its asymptotic value. 
Here 9j is the opening angle of the jet at radius r and 
Fo and 9q are, respectively, the initial Lorentz factor and 
opening angle at the injection radius ro- 

Let us first consider the scenario in which the precur- 
sor is related to the jet breakout (e.g. Ramirez-Ruiz et 
al. 2002a,b; Waxman & Meszaros 2003; Morsony et al. 
2007). In this scenario, the precursor is produced by the 
bow shock emission or the cocoon emission when the jet 
head breaks out of the stellar surface. After this, the 
shocked jet phase emerges, which may have a relatively 
small opening angle for ultra-relativistic material due to 
that the mixing of the shocked jet with stellar material 
lowers its terminal Lorentz factor, so that an observer lo- 
cated at a larger viewing angle would see a dark phase. 
The relativistic reverse shock may be located well below 
the stellar surface at the time when the jet head has just 
reached the surface. One may think that this dark phase 
could last long enough, if the shocked jet material moves 
with a low velocity. However, once the jet head reaches 
the stellar surface, the pressure in front of the jet head 
decreases suddenly and a rarefaction wave will form and 
propagate back into the shocked jet material at the speed 
of sound. The speed of sound in the shocked jet plasma is 
relativistic, c s = c/V3, even if the head velocity Vh <C c. 
When this rarefaction wave arrives at the reverse shock, 
the pressure of the shocked jet material also drops and 
it can no longer decelerate the fast unshocked jet. Sup- 
posing that the reverse shock is located at a position well 
below the stellar surface when the jet head reaches the 
stellar surface, the width of the shocked jet is A <; i?*, 
and the time that the rarefaction wave takes to arrive at 
the reverse shock is t\ — A/c s <; R±/c s — 6R±si s. A new 
system of forward and reverse shocks will form, and the 
rarefied shocked jet plasma is accelerated to a relativistic 
velocity with a Lorentz factor 

r, l2 = r M [4r, ll (i + Cs t/A) 3 ] 1 / 4 (1) 

where T^i is the Lorentz factor of the shocked jet be- 
fore the rarefaction wave propagating back (see Eq.(15) 
of Waxman & Meszaros 2003). Even for a sub-relativistic 
velocity Thi — 1, the rarefied shocked jet is accelerated to 
a relativistic velocity with T^2 J> 2 after a time t\. The 
unshocked jet then takes a time ti ~ R±/c = 3-R^nS 
to reach the stellar surface. Supposing that the internal 
shock occurs at i? 7 where it produces the burst emission, 
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the additional delay relative to the precursor emission is 
£3 = R 1 /2T 2 c. Therefore the total gap period between the 
precursor and the main burst would be only 

T gap =ti+i2+i3-10i?*,iis. (2) 
3. FALLBACK COLLAPSAR SCENARIO 

The above described difficulty in interpreting long (~ 
100 s) precursor gaps through a single jet breakout makes 
it necessary to explore other alternatives. Here we propose 
a model where the gap is related to the central engine ac- 
tivity, and the precursor is produced by an initial weak 
jet launched before the main jet. The <~ 100s long gap is 
rcmisccnt of the natural timescales calculated numerically 
in the fallback collapsar (type II collapsar) model (Mac- 
Fadyen et al. 2001), where the fallback disk forms minutes 
to hours after the initial core collapse 5 . In this scenario, 
the collapse of the iron core initially forms a proto-neutron 
star and launches a supernova shock. However, this shock 
lacks sufficient energy to eject all the matter outside the 
neutron star, especially for more massive helium cores, be- 
cause the gravitational binding energy of the helium core 
increases with mass roughly quadratically, while the explo- 
sion energy decreases (Fryer 1999). If the explosion energy 
is lower than 10 51 erg, as the supernova shock decelerates 
while it travels outward, some of the expanding material 
can decelerate below the escape velocity and falls back. 

During the initial collapse, it is possible that a weak jet 
with an energy of, say, a few times 10 50 erg, is produced as 
a result of MHD processes in the collapsed core, or through 
propeller effects associated with a proto-neutron star, etc. 
(e.g. Wheeler et al, 2000). If such a jet is weak enough 
so as not to immediately disrupt the star, a fallback ac- 
cretion disk can form in the core on a fallback timescale, 
which can launch a stronger jet producing the main burst. 
We argue that such initial weak jets may be responsible 
for precursors with non-thermal spectra, through internal 
dissipation mechanisms such as internal shocks or recon- 
nection, after they exit the stellar progenitor. The non- 
thermal spectra are consistent with those observed from 
most GRB precursors (e.g. Lazzati 2005; Cenko et al. 
2006; Romano et al. 2006; Page et al. 2007). 

After a time ~ 100 s following the initial collapse, in 
the case of a weak supernova, enough gas has fallen back 
for the neutron star to collapse to a black hole. Gas that 
continues to fall back with sufficient angular momentum 
will settle into a disk (MacFadyen et al. 2001). The black- 
hole accretion disk system may then produce stronger jets 
which, through neutrino annihilation or MHD processes 
(e.g. Blandford-Znajek mechanism), can power the main 
burst event. 

4. JET PROPAGATION THROUGH THE STELLAR 
PROGENITOR 

We investigate now whether a weak precursor jet of lu- 
minosity Lj ~ 10 49 ergs _1 lasting for <~ 10 seconds (the 
precursor duration) can break out of the star. Previous 
analytic studies on the jet dynamic inside the star and 
the jet breakout time have been done by, e.g. Meszaros & 
Rees 2001; Ramirez-Ruiz et al. (2002b); Matzner (2003), 
Lazzati & Begelman (2005) and Toma et al. (2006). As 

5 Cheng & Dai (2001) suggested a similar, "two-step" jet model, w 
model of GRBs. 



recently found by Woosley & Zhang (2007) in the simula- 
tion of jet propagation through the star, a lower-energy jet 
takes a longer time to break out of the star. A jet with a 
lower Lorentz factor may also take a longer time, while at 
present there is no constraint yet on the Lorentz factor of 
the precursors jets since their spectra do not display sig- 
nificant ;> MeV emission. Since the precursor lifetime is 
usually only about ten seconds and the breakout time must 
be shorter than this (otherwise the jet will be chocked), it 
is not obvious whether it can really break out of the star 
or not. We will show below that, for a low-energy jet to 
be able to break out of the star in a time shorter than the 
precursor duration, it must have a moderately high asymp- 
totic Lorentz factor, r\ = L/Mc 2 ;> 10. The asymptotic 
Lorentz factor of this weak precursor jet can, however, be 
lower than that of the typical main GRB ejecta, which 
are characterized by r\ > 100, and this could account for 
the softer spectrum of the precursor compared to the main 
burst. Phcnomcnological relations between the energy (or 
luminosity) and the spectrum of the burst have been found 
by e.g. Amati et al. (2002) and Yonetoku et al. (2004), in 
which a lower energy (or luminosity) burst tends to have 
a lower peak energy, hence a softer spectrum. 

4.1. Jet crossing time 

The jet head velocity is given by the longitudinal balance 
between the jet thrust and the ram pressure of material 
ahead of it (Meszaros & Rees 2001; Matzner 2003) 

«*=(a^fe) = 7 x ioM^^V-- 1 (3) 

where 8j is the opening angle of the jet while propagating 
inside the star, whose value will be derived below, and p 
is the stellar density at radius r. The presupernova den- 
sity profile is assumed to be roughly described by p oc r~ 3 
out to the edge of the He core at r = 10 11 cm (MacFadyen 
et al. 2001), at which point we take the density to be 
p = lgcm~ 3 . The total time taken by the jet head to 
move from the interior of the star to the surface is 

= ? = ^LJ^ s (4) 

J Vh Vh 

where a > 1 is an integrating factor over r (discussed af- 
ter Eq.(10)). While the jet is propagating inside the star 
with a sub-relativistic velocity, a significant fraction of jet 
"waste" energy is pumped into the cocoon surrounding the 
advancing jet. The unshocked jet moving with Tj 3> 1 is 
now constrained in the transverse direction by the pres- 
sure of the cocoon, so that it gets collimated and becomes 
more penetrating. The opening angle 9 of the unshocked 
jet is determined by the balance between the pressure of 
the unshocked jet in the transverse direction and that of 
the cocoon. 

The energy stored in the cocoon is roughly E c ~ Lj (t — 
r/c) when the sub-relativistic jet reaches radius r after a 
time t. The volume of the cocoon is V c = | r\r, where 
r± = v±t is the transverse size of the cocoon (Matzner 
2003). Since the cocoon pressure is much larger than that 
of the stellar material, the transverse expansion velocity 

attempts to solve the baryon contamination problem in the hypernova 
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vj_ of the cocoon is given by the balance between the co- 
coon pressure p c and the ram pressure of the stellar mate- 
rial, i.e. v± = (pc/p) 1 / 2 . From p c = E C /3V C , we get 

»-(s?r-(W- («> 

The pressure in the unshocked relativistic jet in the 
transverse direction, from the relativistic Bernoulli equa- 
tion, is pj t ± = ^Lj/(2Ttr 2 6 2 cT 2 ). As the internal energy 
gets converted into the kinetic energy of the jet, Tj in- 
creases. Depending on whether the jet Lorentz factor has 
saturated or not before the jet breaks out of the star, there 
are two cases for the evolution of I^-: one is the high rj 
case (case I), in which Tj oc (r9j) = To(r /ro)(9j /6q) < rj 
all along the path (see the left panel of Fig. 1) , and the 
other is the low 17 case (case II) in which Tj has reached 
its saturation value rj at some radius inside the star (see 
the right panel of Fig. 1). The latter case is more easily 
satisfied when 77 is lower or when the initial injection ra- 
dius ro is smaller. The jet pressure is different in the two 
cases, being given by, respectively, 

f Ljrlel/^cr^Tl) Tj < rj (Case I); 
ft '- L \ L 3 /(6Trcr 2 6 2 rj 2 ) Tj = rj (Case II) { ' 

From pj t ± = p c , we get the opening angle for the two cases 
respectively, 

„ _ f ( a 2 V6487rpc 3 ) 1 /i4( ro Vr )4/7 r -5/7. m 
3 ~ \ (a 2 L j /6487rr 2 pc 3 ) 1 / 6 ??- 4 / 3 ( ' 

For case I, with the parameters #0 = 1, r = 1, r = 10 7 cm 
and the stellar density profile p = l(r/10 n cm)~ 3 gem -3 , 
the opening angle is 

6j = 1.5 x lO-V/^r^rS?, (8) 
while for case II, 

^^SxlO-V/^r^V 473 - (9) 

Note that these values of the opening angles are the ones 
at the jet breakout time, which are not the final values that 
we expect in the internal shock phase or in the afterglow. 
After exiting the star, the hot jet is expected to expand 
freely to an opening angle comparable to 1/Tb r ~ 0.1, 
where Tb r is the Lorentz factor of the jet at the breakout 
time. 

With the above values of the opening angles, we can ob- 
tain the jet breakout time from Eq.(4) for the two cases, 

_ f a^r^C-^)^(^/7 = 2a^LT% 7 r 4 / !7 7 S ; 

^ " { a 4/3 r 5/3 ( ^ ) -l/3 ?? -4/3 = 4^/3^/3^-1/3^ 4/3 s 

(10) 

Here a is the integrating factor over r in Eq.(4). In case 
I, since p oc r~ 3 starts roughly above r = 10 9 cm, be- 
low which the free fall density profile p ex r -3 / 2 may ap- 
ply (Meszaros & Rees 2001), a =1 ln(r/10 9 cm) ~ 4 for 
r = 10 11 cm. In case II, the integration over r in Eq.(4) 
gives a 2 = 3/2. We can see that a low-luminosity, weak jet 



takes a longer time to break out of the star, consistent with 
recent findings by Woosley & Zhang (2007) from the nu- 
merical hydrodynamic simulation of the jet propagation. 
The dependence of the breakout time ibr on the luminosity 
is illustrated in the inset plot of the left panel of Fig. 1. Re- 
quiring that the breakout time ibr is less than the duration 
of the precursor T p ~ 10 s, we get, for case I, 

Lj * 10 49 (a 1 /4) 8 / 3 r 4 / 7 3 T- 1 7 / 3 ergs- 1 . (11) 

The asymptotic Lorentz factor rj in this case is required 
to be larger than Tb r = To(rdj /Vo#o) — 15, from the high 
Tj assumption of case I. For case II, the breakout time is 
sensitive to rj, more than to the jet luminosity, and we get 

,,^10(a 2 /1.5)r 1 1 / 2 L7V 4 T- 1 3 / 4 (12) 

by requiring the breakout time to be shorter than T p ~ 10s 
(see the inset plot of the right panel of Fig.l). Therefore 
we see that in both cases, it is only when rj ;> 10 that the 
jet breakout time can be <J 10 seconds, for a luminosity of 
the order of Lj ~ 10 49 ergs _1 . This reflects the fact that 
a lower Lorentz factor jet is less penetrating and needs a 
longer time to break out of the stellar progenitor. This 
shows that even the precursor jet should have a relatively 
high Lorentz factor, although it need not be as large as 
;> 100 as in the case of the main burst. 

4.2. Can the precursor channel close up again? 

After the weak jet breaks out, the cocoon material es- 
capes. There is now a wide funnel opened by the weak 
jet and its accompanying cocoon along its path. As evi- 
denced by GRB050820A (Cenko et al. 2006), a burst with 
precursor, the main burst jet is collimated to 9 ~ 0.1, 
like other typical long bursts without precursors. Since a 
closed channel is a necessary condition for the main jet 
to be collimated, one may ask whether the initial chan- 
nel, made by the precursor jet, will have had time to close 
again by the time that the main jet starts to make its way 
out, ^100 seconds later. 

As the weak jet propagates outward, the cocoon pushes 
the stellar material sideways. The opening angle of the 
cocoon, which is also the funnel opening angle after the 
cocoon flows out, is approximately given by 

»/-«•- T-(^r-«^'- < 13 > 

where the last equality is obtained with the help of the 
second line of Eq.(10). For Tb r — 10, the cocoon opening 
angle is about tens of degrees. 

The cocoon can drive a transverse shock into the stellar 
material along its passage, heating it up, so the pressure 
of the shocked stellar material (p e ~ pv\) is comparable 
to that of the cocoon, i.e. p e ~ p c . After the cocoon flows 
out, this external stellar material will re-expand back into 
the funnel at a speed of sound given by 

c s = VPe/P - VPc/p = v±. (14) 

This means that in a time comparable to the precursor jet 
breakout time, the channel closes up again. From Fig. 5 of 
Lazzati (2005), it is known that the precursor duration is 
always shorter than the gap period of time. Therefore, we 
conclude that the channel opened by the precursor jet is 
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indeed closed again before the main jet starts to propagate 
out. 

The free-fall time of the material along the rotation axis 

is 

-1/2 



t g = 3000 



1 gem 



s. 



(15) 



For the C/O core material of the collapsed progenitor, 
whose density is p <~ 500gcm~ 3 , its free-fall time is compa- 
rable to the ~ 100 s gap period and therefore this material 
will be essentially evacuated. However, for He core mate- 
rial with p ~ 1 gcm~ 3 , the free-fall time is much longer 
than the gap period, and the funnel in this part (aside 
from reclosing) remains almost the same as before the col- 
lapse, by the time that the main jet starts to propagate 
out. 

4.3. Transient thermal emission from the open funnel 

The channel made by the initial precursor jet can, while 
it remains open, produce an interesting emission signature. 
As the cocoon expands transversally into the stellar mate- 
rial, it will drive a shock. Since this shock takes place in 
the highly optically thick interior of the progenitor star, it 
will become a radiation-dominated shock, similarly to the 
case of the supernova shock propagating inside the star 
(Weaver 1976). The velocity of the shock driven by the 
cocoon into the He core material can be derived from v s = 



(Pc/P) 



1 / 2 = vj_, where p c is given by Eq.(5). This gives a 
shock velocity v s ~ 10 9 cms~ 1 L^4g^y i 4 r^ 1 3 ^ 4 p _3 / 8 a~ 1 / 4 . 
For a radiation-dominated shock, this corresponds to a ra- 
diation temperature of 



Ljpv h 



1/8 /%\V 



'r n 1/4 KeV 



(16) 

where a is the Boltzmann energy density constant of radi- 
ation. 

To check the self-consistency of having assumed a 
radiation-dominated shock, we need to compare the the 
radiation diffusion length A 7 = c/ (3n e v s a^) (where n e 
is the electron number density of the plasma) with the 
stopping distance of the ions due to Coulomb collisios I. 
According to Weaver & Chapline (1974), the Coulomb 
friction is the dominant dissipation mechanism control- 
ling the shock width (i.e. the shock is collisional), pro- 
vided that the stopping distance of the ions I is larger 
than the radiation diffusion length A 7 . Otherwise, the 
shock is radiation-dominated. For an electron tempera- 
ture T e = T r , the ion-electron collision rate is Vi e = 2 x 
10 10 (p/lgcm- 3 )(T e /10KeV)- 3 / 2 s- 1 , while the ion-ion col- 
lision rate, for an ion temperature of T, = (3/16)m p ?; 2 = 
200KeV, is v u = 1 x 10 10 (p/lgc m - 3 )(T i /200KeV)- 3 / 2 s- 1 
(Spitzer 1962; Waxman & Loeb 2001). Thus, the stop- 
ping of the ions is dominated by interaction with the 



radiation temperature of T r ~ lOKeV, Compton scatter- 
ing gives a cooling rate of t^l mp — (8crr/3m e c)aT 4 = 
8 x 10 10 (T,./10KeV) 4 s- 1 , while the Bremsstrahlung cool- 
ing rate is t'^ rem = 3x 10 8 (p/lgc m - 3 )(T e /10KcV)- 1 / 2 s- 1 . 
Thus, Compton scattering also provides the dominant 
mechanism for cooling the electrons, and this can re- 
duce the ion temperature significantly from those that 
would occur in a Coulomb viscosity-dominated shock. 
These photons diffuse ahead, creating sufficient pressure 
to decelerate the electrons (as viewed in the frame of the 
shock front) and the ions are then decelerated by these 
electrons through Coulomb friction. We note that the 
above rough quantitative argument is only an approxi- 
mate treatment of this issue. A more detailed numerical 
treatment taking into account the energy and momen- 
tum balance in the shock transition as well as relativistic 
effects 6 , as performed by Weaver (1976), is beyond the 
scope of the present work. However, there is a useful 
numerical criterion for checking the self-consistency of a 
radiation-dominated shock. According to Weaver (1976), 
the radiative heat transport suffices for producing the nec- 
essary dissipation in a strong shock wave if the ratio of the 
radiation pressure to matter pressure P r /P m exceeds 4.45. 
This condition is satisfied in our case, the ratio of the two 
pressures being P r /P m = aT^/n e kT e = pv 2 /n e kT r ~ 100 
for v s — 10 9 cms _1 and T e = T r = lOKeV, justifying the 
assumption of a radiation-dominated shock. 

The shock-heated walls of the jet channel, before the lat- 
ter closes up again, can radiate a transient thermal pulse 
with characteristic photon energy kT ~ lOKeV and a lu- 
minosity 



L x ~ o-T 4 ir(R+6 f ) 2 

= 2 x 10 48 a 3 /%V 4 2 9 (^Xn^ergs- 1 



(17) 



electrons, with a mean free path given by U e = v ie v s = 
0.05cm(p/lgcm- 3 )- 1 (T e /10KeV) 3 / 2 (i; iS /10 9 cms- 1 ). The 

radiation diffusion length is A 7 = 20cm(p/lgcm _3 ) _1 (i; s /10 9 cmBlo^ni.2006; Metzger et al. 2007). Following Wheeler et 
Since A 7 3> I 



As the cocoon escapes from the channel, the cocoon 
pressure decreases significantly, as does the radiation en- 
ergy density of the photon field. Finally, as the pressure 
wanes, the channel closes up again under the pressure of 
the external stellar material. The emission from the 
channel walls vanishes once the channel has reclosed, so 
this transient thermal pulse is expected to last for a time 
comparable to the jet crossing time given by Eq.(10). This 
could appear as a thermal precursor of ~ 10 seconds du- 
ration, besides the non-thermal emission component from 
the cocoon itself. 

5. POSSIBLE FORMATION MECHANISMS OF THE 
PRECURSOR JETS 

How does the weak precursor jet form during the initial 
core collapse? There is a long-standing speculation that 
when rotation and magnetic fields are taken into account, 
the core collapse of massive stars can lead to some form of 
MHD outflows, although it is unknown whether they can 
power supcrnovac or gamma-ray bursts (e.g. LeBlance 
& Wilson 1970; Wheeler et al. 2000, 2002; Woosley & 



the scale (width) of the shock is deter- 
mined by A 7 . Except for a thin layer at the leading edge 
of the shock front, the Coulomb friction is unimportant, 
and the shock is effectively radiation-dominated. For a 



al. (2000), here we consider that a relatively weak jet of 
energy a few times 10 50 erg is produced by a rotating proto- 
neutron star during the initial core collapse phase. This is 
motivated by the consideration that the rotation energy of 



6 For a radiation temperature T r ~ lOKeV, the pair production effect is unimportant (Weaver 1974) 
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a proto-neutron star with a radius of 50 km is suitable for 
GRB precursor jets. Immediately after the core collapse, a 
proto-neutron star (PNS) forms with a radius of <~ 50km 
and mass ~ 1.4Mq. For a specific angular momentum 
of the progenitor core of j ~ 10 cm s , the angular 
rotation velocity of the PNS is ft ~ 400s _1 , assuming an- 
gular momentum conservation. The rotation energy of the 
proto-neutron star is 

(18) 

Over a cooling time ~ 5 — 10s, the hot proto-neutron star 
dclcptonizes by neutrino emission and contracts to form 
the final neutron star. If a fraction, e.g. <~ 10% of the 
rotation energy is tapped by some MHD process and con- 
verted into a relativistic outflow, this outflow will have 
the right energy and timescale for a GRB precursor. One 
possible scenario for such an MHD process is that the dif- 
ferentially rotating neutron star can wind up the poloidal 
seed field inside the proto-neutron star into strong toroidal 
fields, which then emerge from the star through buoyancy 
(Kluznizk & Ruderman 1998; Dai & Lu 1998; Dai et al. 
2006). Once the field emerges through the stellar sur- 
face, it will trigger an explosive reconnection and release 
the magnetic energy, resulting in a high-entropy fireball. 
Assuming a critical field value required for buoyancy of 
Bf ~ 10 15 G corresponding to the density of the proto- 
neutron star (Wheeler ct al. 2002), the number of revolu- 
tions required is nj = Bj/2ttBo = 160Bj^ 5 Bq1 2 where B 
is the magnetic field of the proto-neutron star. The ampli- 
fication timescale before the field is expelled by buoyancy 
is thus Tf ~ %-Pp N s — 2.5B/ i i5B^"i 2 s, where Ppns IS the 
spin period of the proto-neutron star. 

Note that the infalling matter from the exterior of the 
proto-neutron star is insufficient to suppress this buoyant 
magnetic field pressure in a proto-neutron star. The ram 
pressure at i?pNS = 50km for a mass infall rate M is 

(19) 

(Fryer 1999), which is much smaller than the magnetic 
field pressure ps = Bj/8ir ~ 4 x 10 28 Bj 15 ergs _1 , where 

Vff = (2GM/i?PNs) 1/2 is the free-fall velocity , M is the 
mass infalling rate in units of Mqs -1 and M is the proto- 
neutron star mass. The infall rate of matter at early times 
(the first 10 s) in the fallback collapsar is uncertain, but 
as a rough estimate, the accretion rate (~ O.OIMqs -1 ) 
during the fallback phase (the first 100 s) can be used (see 
Fig. 5 of MacFadyen et al. 2001). However, after about 
10 s neutrino cooling will tend to shrink the proto-neutron 
star size, and since the ram pressure is proportional to 
i? -5 / 2 , even if a magnetar with B <^ 10 15 G could form af- 
ter contraction, its activity would be suppressed, due to its 
smaller size R, by the ram pressure of the infalling matter 
after the infall time <~ 100 s. 

Since the activity of a central neutron star with a mag- 
netic field B <, 10 15 G can be suppressed by the infalling 
matter, no significant emission would be seen after the 
proto-neutron star phase, until the accretion disk forming 
at a later time results in the formation of a newborn black 



hole. In the fallback collapsar scenario, the accretion rate 
depends on the initial supernova shock energy. For a weak 
SN explosion with an energy <; 0.5 x 10 51 erg, the accre- 
tion rate is about <~ O.OIMqs -1 (see Fig. 5 of MacFadyen 
et al. (2001)). With the above accretion rate, the neutron 
star will collapse to a black hole on a timescale of the or- 
der of 30 — 100 seconds. After the black hole forms, gas 
that continues to fall in with sufficient angular momentum 
will settle into a disk, which can power a main jet through 
MHD processes or neutrino annihilation in the disk-black 
hole system. This ~ 100 s accumulation timescale of the 
accretion material may be responsible for the long time 
delay of inactivity seen in some bursts with precursors. 

Another possible mechanism for causing a long-time gap 
might be the development of instabilities in the hyperac- 
cretion disk (e.g. Perna et al. 2006; Janiuk et al. 2007). 
Janiouk et al. (2007) found that the opaque inner re- 
gions of the disk become unstable to convective motions. 
This may give rise to the short-term time variability seen 
in the prompt phase, but can not explain the long-time 
gap. The long time delay may require, e.g., that the disk 
is gravitationally unstable, fragmenting into rings. The 
gravitational instability sets in when the Toomre param- 
eter (Toomre 1964) Q — c s k/ttGT, < 1, where c s is the 
sound speed, k is the local epicyclic frequency and £ is 
the mass surface density. The disk in the case of a merger 
of compact objects is found to be gravitationally stable, 
since Q > 1, according to Lee & Ramirez- Ruiz (2002) and 
Lee et al. (2005). In order to explain the X-ray flares 
seen in short GRBs, Perna et al. (2006) speculated that 
the outer part of the disk may become gravitationally un- 
stable, when the accretion rate is as high as 10M Q s _1 . 
However, for fallback collapsars, the outer radius at which 
this criteria is satisfied may lie well outside the the ac- 
cretion shock radius seen in the fallback calculations by 
MacFadyen et al.(2001), which casts into doubt the via- 
bility of this mechanism for the long-time gap. 

The long gap may also be produced by a modulated 
relativistic wind without postulating any quiet phase in 
the central engine (Ramirez- Ruiz et al. 2001). It was also 
pointed out that an approach to test this would be through 
observations of the prompt afterglow emission from the re- 
verse shock. 

6. DISCUSSION 

We have explored the constraints posed by GRB pre- 
cursors with long time gaps on the jet dynamics while 
propagating through the progenitor stars and jet param- 
eters. We showed that a ~ 100 s long gap is hard to 
interpret in terms of a single jet breaking out through the 
stellar envelope, and we suggested a possible two stage 
scenario, involving a precursor jet powered by an initial 
proto-neutron star, and a main burst powered by a fall- 
back collapsar leading to a black hole. In this two stage 
central engine scenario, the precursor jet in the first stage 
must be weak, so as not to disrupt the star, otherwise 
there would be no fallback accretion disk to power the 
second stage main jet. Since the observed precursor du- 
rations are usually short, of the order of ten seconds, the 
jet propagation time through the star poses severe con- 
straints on the jet physical parameters, such as the jet 
luminosity and Lorentz factor. One of our aims has been 
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to give an analytic calculation of the jet propagation time, 
in order to explore the properties of precursor gap time 
models. Comparison of our analytical results with numer- 
ical hydrodynamical simulations such as those of Woosley 
& Zhang (2007) show that the analytical dependence of 
the propagation time on the jet luminosity is in qualita- 
tive agreement with these simulation results. 

We have studied the effects of a precursor jet breaking 
out from the progenitor and its related possible observa- 
tional manifestations, within the context of our two-stage 
scenario. The precursor jet opens a channel through the 
star as it makes its way out. The channel will re-close 
after the precursor jet and its cocoon have escaped, on a 
timescale of ~ 10 s. This allows a main jet to form as a 
result of a collapsar black hole accretion, if the precursor 
is weak enough so as not to disrupt the star, the main jet 
being collimated by the stellar pressure and propagating 
outward in the usual manner. These analytical results are 
approximate, and numerical hydrodynamical simulation 
studies of such a two-stage jet model would be required 
to verify its behavior in more detail. 

We have considered some of the possible formation 
mechanisms for the precursor jet in the framework of the 
two-stage fallback collapsar scenario, involving MHD ac- 
tivity of an initial proto-neutron star. We stress that a 
proto-neutron star MHD outflow is just one possibility, es- 
pecially considering the many uncertainties involved both 
in the proto-neutron star and in the fallback collapsar sce- 
narios. The possible origin of the ~ 100 s inactivity of the 
central engine, in this two-stage fallback collapsar frame- 
work, is also tentative. However, for any possible two-stage 
scenario, the above discussion on the jet-star interaction 



dynamics remains applicable, and the constraints we have 
discussed remain applicable. 

In summary, we have argued that the long time (~ 100 
s) gap between the precursor and the main burst may be 
attributed to the central engine behavior. One natural 
scenario resulting in a ~ 100 s delay is the fallback collap- 
sar model (MacFadyen ct al. 2001). We suggest that the 
initial precursor jet may arise in an initial stage of the col- 
lapsed core, possibly related to MHD processes or neutrino 
emission occurring in a proto-neutron star. We have calcu- 
lated the crossing time of the weak precursor jet and have 
shown that this initial weak jet can break out through the 
progenitor envelope in a time less than ten seconds, which 
is a typical precursor duration, provided that it has a rela- 
tively high Lorcntz factor T J> 10. The channel opened by 
the precursor jet along the rotation axis of the collapsar 
may produce a transient weak thermal precursor, before it 
recloses due to the re-expanding stellar material which was 
heated by the cocoon of the precursor jet. This provides 
the conditions required for the collimation of the main jet, 
which starts out after a typical fall-back timescale of ~ 100 
seconds, producing a black hole leading to the main GRB. 
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NAG5-13286, NSF AST 0307376, and the National Natu- 
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National Excellent Doctoral Dissertations of China (for 
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Fig. 1. — The evolution of the jet Lorentz factor with radius for the two different cases defined in the text. The left panel is the high 
entropy case, where the jet Lorentz factor is not saturated before breaking out of the star and the right one is the lower entropy case, where 
the jet Lorentz factor has already saturated before breaking out. The inset figures in the two panels show the dependence of the jet breakout 
time t br on other parameters for the two cases respectively. Here ro is the initial injection radius, To is the initial injection Lorentz factor 
and R t is the stellar radius of the progenitor. 



